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Corrections of orders a 5 and a 6 are calculated in the hyperfine splittings of IS 
and 2S - energy levels in the ion of muonic helium. The electron vacuum polariza- 
tion effects, the nuclear structure corrections and recoil corrections are taken into 
account. The obtained numerical values of the hyperfine splittings -1334.56 meV 
(15 state), -166.62 meV (2S state) can be considered as a reliable estimate for 
the comparison with the future experimental data. The hyperfine splitting interval 
A12 = {8AE h f s (2S)-AE h f s (lS)) = 1.64 meV can be used for the check of quantum 
electrodynamics . 
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I. INTRODUCTION 

The ion of muonic helium (n\He) + is the bound state of the negative muon and helion 
(iHe). The lifetime of this simple atom is determined by the muon decay in a time r M = 
2.19703(4) x 1CT 6 s. An increase of the lepton mass in muonic hydrogenic atoms as compared 
with electronic hydrogenic atoms (m^/rrie = 206.7682838(54) |l|) leads to the growth of three 
effects in the energy spectrum: the electron vacuum polarization, the nuclear structure and 
polarizability, the nuclear recoil effect. The first effect is important for the spectrum of 
muonic helium ion ([i\He) + because the ratio of the Compton wavelength of the electron 
to the Bohr radius of the atom /iZa/m e « 1.45415 is sufficiently close to the unity. The 
second effect of the nuclear structure is of the utmost importance because the muonic wave 
function has a significant overlap with the nucleus. Finally, the increase of the recoil effects 
is determined by the ratio of the muon and helion masses m^/m(^He) ~ 0.0376 [l|. Despite 
the fact that this ratio is small, it is many times larger than the fine structure constant a 
(a -1 = 137.03599911(46) [l|). Moreover, some recoil effects contain the peculiar logarithm 
of the muon to helion mass ratio ln(m(|.ffe)/m /Lt ) « 3.28 what leads to the numerical growth 
of the contribution. 

High sensitivity of the bound muon characteristics to the distributions of the nuclear 
charge and magnetic moment in light muonic atoms (muonic hydrogen, ions of muonic he- 
lium) can be used for more precise determination of the charge radii of the proton, deuteron, 
helion and a - particle @, S, 0] ■ Moreover, the measurement of the hyperfine structure of 
light muonic atoms allows to obtain more precise values of the Zemach radii and to improve 
the accuracy of the theoretical calculations of the hyperfine splittings in the corresponding 
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electronic atoms. 

Theoretical investigations of the energy levels in light muonic atoms (the nuclear charge 
Z = 1,2) were carried out many years ago in Refs. [!, @, 0, @, 0] on the basis of the Dirac 
equation. In the calculations of the energy intervals (2P 3 / 2 — 25i/ 2 ), (2Pi/ 2 — 25*1/2) in the 
ions of muonic helium (|Pe, \He) the different type corrections were considered with the 
precision 0.01 meV. The energy transitions (25* — 2P) in the ion of muonic helium (y^He) 
were calculated in [1(3, ll[ with regard to the hyperfine structure on the basis of the Dirac 
equation with the accuracy 0.1 meV. 

Although muonic atoms (fi\H), (/j^H), (/i 2 Pe) + , (/i 2 ife) + could be used for yet another 
check of quantum electrodynamics, the experimental study of the energy levels of these 
atoms lags behind the theory. It is appropriate at this point to recall the experiment for the 
Lamb shift (2P — 2S) measurement in muonic hydrogen which is carried out many years 
at PSI (Paul Sherrer Institute) [H, Hf but didn't give while the positive result with the 
necessary accuracy 30 ppm. The only experiment with the successful end was performed on 
the muon beam at CERN 0, EEl with muonic helium (/i 2 Pe) + . Two resonance transitions 
with the wavelengths 811.68(15) nm and 897.6(3) nm corresponding to the fine structure 
intervals (2P 3 / 2 — 2S'i/ 2 ) and (2P^ 2 — 2Si/ 2 ) were observed. In a later experiment [16J the 



resonance transition in the range 811.4 < A < 812.0 nm was not observed. So, at present 
new measurement of the Lamb shift and the hyperfme structure in the atoms (/i 2 Pe) + , 
(/i 2 Pe) + is required. 

There is a need to remark that in the last years the accuracy of the theoretical inves- 
tigations of the energy spectra of simple atoms was increased essentially [sj. New QED 
corrections of order a 6 and a 7 in the energy spectra of muonium, positronium, hydrogen 
atom and ions of electronic helium were calculated pj. For a number of hydrogenic atoms 
(hydrogen atom, ions of helium) the comparison of results of the theoretical investigations 
in QED with the experiment is difficult because the theoretical error in the calculation of 
the nuclear structure and polarizability contributions both for the Lamb shift and hyperfine 
structure remains very large and exceeds considerably the experimental errors. The progress 
in this field can be achieved due to new experimental investigations of the structure and 
polarizability of the proton and other nucleus and the use of light muonic atoms. 

It is important to keep in mind that all contributions to the energy spectra of light muonic 
atoms can be divided into two groups. The corrections of the first group were obtained in 
the analytical form in the study of the energy levels of muonium, positronium and hydrogen 
atom. The second group includes numerous corrections of the electron vacuum polarization, 
the nuclear structure, recoil effects which are specific for each muonic atom. The aim of 
this work consists in the analytical and numerical calculation of corrections of orders a 5 
and a 6 in the hyperfine structure of S-states in the muonic helium ion (/i 2 Pe) + on the 



basis of the quasipotential method in quantum electrodynamics [17|, uM- We consider such 
effects of the electron vacuum polarization, recoil and nuclear structure corrections which 
have the crucial importance to attain the high accuracy of the calculation. Numerical values 
of corrections are obtained with the precision 0.001 meV. So, the purpose of our investi- 
gation consists in the improvement of the earlier performed calculations 0, Q, H, 0| of the 
hyperfine splitting in the ion of muonic helium and derivation of the reliable estimates in the 
hyperfine structure of 15 and 2S'-states which could be used in conducting a corresponding 
experiments. Modern numerical values of fundamental physical constants are taken from 
the paper [l[: the electron mass m e = 0.510998918(44) ■ 10~ 3 GeV, the muon mass m M = 
0.1056583692(94) GeV, the fine structure constant a' 1 = 137.03599911(46), the helium mass 
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m{lHe) = 2.80839142(24) GeV, the helium magnetic moment fi h = -2.127497723(25) in 
the nuclear magnetons, the muon anomalous magnetic moment a M = 1.16591981(62) • 1CT 3 . 

II. EFFECTS OF ONE-LOOP AND TWO-LOOP VACUUM POLARIZATION IN 

THE ONE-PHOTON INTERACTION 

Our approach to the investigation of the hyperfme structure (HFS) in the muonic helium 
ion is based on the quasipotential method in quantum electrodynamics [H, 20, 21], where the 
two-particle bound state is described by the Schroedinger equation. The basic contribution 
to the interaction operator of the muon and helion for ^-states is determined by the Breit 
Hamiltonian (22|: 

H B = H + AV£ S + AV* S \ H = f-—, 
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where m 1; m 2 are the muon and helion masses, m p is the proton mass, fih is the helion 
magnetic moment. The potential of spin-spin interaction (3) gives the main contribution 
to the energy of hyperfme splitting of 5-states (the Fermi energy). Averaging (3) over the 
Coulomb wave functions of IS and 2S states 
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we obtain the following result (the difference of the triplet and singlet states): 



AE h / s (nS) = 8 ^ Z3a \ h 



IS: - 1370.725 meV 



3mim p n 3 \ 2S : — 171.341 meV 



(4) 
(5) 

(6) 



The muon anomalous magnetic moment does not enter in Eq.(6). The correction of 
the muon anomalous magnetic moment (AMM) to the hyperfme splitting is conveniently 
represented by the use the experimental value a M = 1.16591981(62) • 10~ 3 [if: 



AE h J u s (nS)=a,AE h / s (nS) 



IS : - 1.598 meV 
2S : - 0.200 meV' 



(7) 



The contribution of the relativistic effects of order a 6 to the hyperfme structure is also 
known in the analytical form (3[: 
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One-loop electron vacuum polarization correction in the interaction operator is deter- 
mined by the following expression in the coordinate representation 22] : 



8anh CjfTj a 
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where = y/t; 2 — 1(2£ 2 + l)/£ 4 . To obtain (9) it is necessary to use the following replace- 
ment in the photon propagator: 
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Averaging (9) over the wave functions (4) and (5), we find the correction of order a 5 to the 
hyperfine splitting: 
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-4.203 meV, 

(11) 
(12) 
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Changing the electron mass m e to the muon mass nil in Eqs. (11), (12), we obtain the 
muon vacuum polarization contribution to the hyperfine splitting of order a 6 because the 
ratio W/m>i <C 1. The corresponding numerical values are included in Table 1. The contri- 
bution of the same order a 6 is given by the two-loop electron vacuum polarization diagrams 
in Fig.l(b,c,d). To derive the interaction operator corresponding the amplitude with two 
sequential loops in Fig. 1(b) we use the double change (10). In the coordinate representation 
the interaction operator has the form: 
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Corresponding correction to the hyperfine splitting of IS" and 2S levels can be presented 
as a three-dimensional integral over variables r, £ and 7]. After that the integral over r is 
calculated analytically and over £, 77 numerically with the result: 
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In a similar way we calculate the two-loop vacuum polarization contribution of order a 6 
shown in Fig.l(c,d). In this case the potential is determined by the relation: 
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where the function 

f(v) = v{(3-v 2 )(l+v 2 ) 
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Lii{z) is the Euler dilogarithm. Numerical contributions of the operator (16) to the HFS 
are included directly in Table 1. The role of the vacuum polarization effects in the HFS of 
the muonic helium ion extends further. There exists a number of contributions where the 
electron vacuum polarization enters the potential together with the nuclear structure, recoil 
and relativistic effects in the second order perturbation theory. 



abed 
FIG. 1: Effects of one-loop and two-loop vacuum polarization in the one-photon interaction. 



III. EFFECTS OF ONE-LOOP AND TWO-LOOP VACUUM POLARIZATION 
IN THE SECOND ORDER PERTURBATION THEORY 

The second order perturbation theory (PT) corrections in the energy spectrum of 
hydrogen-like system are determined by the reduced Coulomb Green function G [23j, whose 
partial expansion has the form: 

G n (r,r>) =Y J 9m{ry)Y l MYL^)- (18) 

l,m 
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The radial function g n i(r,r') was presented in Ref. 23j in the form of the Sturm expansion 
in the Laguerre polinomials. The basic contribution of the electron vacuum polarization to 
HFS in the second order PT can be presented as follows (see Fig. 2 (a)): 



A^sopt vp i = 2 < ■ G ■ AV^ >, 

where the modified Coulomb potential 
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The potential AV^ /s (r) is proportional to S(r). So, we need the reduced Coulomb Green 
function with one zero argument. It was derived by means of the Hostler representation as 
a result of the subtraction of the pole term in Refs. (2I. 25]: 
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where C = 0.5772... is the Euler constant, x = Wr. Using Eqs.(21), (22) we can present 
necessary corrections in the HFS of the ion (fi 2 He) + in the form: 

AE&pAlS) = -A££ /s (l,S)^(l + [°° p(£)d£ [°° e- 2 < 1+ ^) g ls (x)dx = -9.260 meV, 

Sir Ji Jo 

(23) 



AE^^S) = AE h / s (2S)^(l + a») f°° p(t)d£ He 

3n Ji Jo 



■x(l- 



'2 rn , g 



;l J 92s(x)dx = —0.869 meV. 

(24) 

The factor (1 + a M ) is introduced in Eqs.(23), (24) so that these expressions contain correc- 
tions of orders a 5 and a 6 . 

Two-loop contributions in Fig.2(b,c,d,e) are of order a 6 . Let us consider first contribu- 
tion which is determined by the potentials (9), (20), the reduced Coulomb Green functions 
(21), (22) and the reduced Coulomb Green functions with both nonzero arguments. It is 



convenient to use the compact representation for it which was obtained in Refs. [24|, |26 
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FIG. 2: Effects of one-loop and two-loop vacuum polarization in the second order perturbation 
theory (SOPT). The dashed line represents the Coulomb photon. G is the reduced Coulomb Green 
function. 



Substituting Eqs.(9), (20), (25) and (26) in Eq.(19), we obtain two contributions for each 
energy level IS* and 2S: 
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While contributions (27), (28) and (29), (30) are individually divergent, but their sum is 
finite. Corresponding numerical values are represented in Table 1. Corrections from two 
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other diagrams in the hyperfine structure can be calculated by the relations (23) and (24), in 
which the replacement of the potential (20) by the following potentials should be performed 

AV v C p _ VP (r) = (^) 2 jf f PWV (-^) ^ (eV^-^ 2 - ) , 

(31) 

ATrr . , 2Za 3 f 1 f(v)dv -^p£= , , 

Omitting further intermediate expressions, which have the general structure quite similar to 
(23), (24), we include in Table 1 numerical values of corrections from potentials (31), (32). 



IV. NUCLEAR STRUCTURE AND RECOIL EFFECTS 



The main contribution of the nuclear structure to HFS of the S'-energy levels including the 
Zemach correction, is determined by two-photon exchange diagrams (see Fig.3). We consider 
that the nuclear charge and magnetic moment are distributed in the space with definite 
densities. The vertex operator of the nucleus 2 He contains the electric Ge and magnetic 
Gm form factors which determine the interaction of the nucleus with the electromagnetic 
field. To calculate the nuclear structure correction of order a 5 we use the equation from 
Ref.(2oT] (the abbreviation "str" is used to designate the nuclear structure correction): 
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To remove infrared divergence in Eq.(33) we must take into account the contribution of the 
iterative term of the quasipotential to the HFS of the atom (n\He) + : 
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where the angle brackets denote the averaging of the interaction operator over the Coulomb 
wave functions and the index "hfs" is indicative of the hyperfine part in the iterative term 
of the quasipotential. Vi 7 is the quasipotential of the one-photon interaction, G' is the free 
two-particle propagator. The integration in Eqs.(33) and (34) can be done by means of the 



dipole parameterization for the Pauli form factor F\ and the Dirac form factor F 2 [271 . |28 
The parameter for this parameterization A 2 can be related with the nuclear charge radius 
r^r: A 2 = 12/r 2 N . Numerical value of rjy = 1.844 ± 0.045 fm is taken from Ref. [8J]. Since 
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the dependence on the principal quantum number n in Eq.(33) is determined by the factor 
1/n 3 , the numerical values of the nuclear structure corrections for the states IS and 2S* 

hfs _ (IS: 48.376 meV 
AEstr -\2S: 6.047 meV (35) 

are cancelled in the special hyperfme splitting interval A 12 = (8AEg/ r s (2S) — AEg/ r s (lS)). 
So, the calculation of the interval A 12 is free of the uncertainty connected with the nuclear 
structure at least in the leading order. The value of the correction (33) is dependent on the 
form of distributions Ge,m- The replacement of the dipole parameterization by the Gaussian 
model leads to the change of the numerical value (35) by 2 per cent. 



a b 

FIG. 3: Effects of the nuclear structure of order a 5 . 

The sixth order over a contribution to the HFS shown in Fig. 4 contains both the nuclear 
structure and the vacuum polarization effects. Using the substitution (10) and Eq.(33) it 
can be written as follows: 

lhfs 2a(Za)V f%/ nu, f 1 v " i 1 ~ i) dv 



AE r s l t r ; vP = ^ '-f- / V VP (k)dk / *' . 2 , (36) 

str ' v ^ m 1 m 2 -K 2 n 3 Jo Jo k 2 (l - v 2 ) + Am 2 e 

where the potential V VP (k) differs on the expression V(k) in Eq.(33) by the additional factor 
k 2 . The amplitude contribution (36) to the energy spectrum should be augmented by two 
iterative terms which are presented in Fig.4(c,d): 

A£fter,str VP = -2<V C xG f x AVyp > hfs = -2 < V$ P x G f x AVp fs > hfs = (37) 



F m e 7r 2 Jo Jo k 2 (l — v 2 ) + 1 
Numerically the sum of corrections (36) and (37) is equal 

a i? h f s a i? h f s \ ^ : 0-^60 meV 

AE stryP + 2AE tMr vp = { 2S: Q Q95 meV (38) 

There exists another contribution of the nuclear structure in the second order PT, which 
is determined by the hyperfme term of the Breit Hamiltonian and the operator of the one- 
photon interaction (see Fig. 5(b)) 

A4 = ^(r). (39) 
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FIG. 4: Nuclear structure and vacuum polarization effects of order a 6 . The dashed line represents 
the Coulomb photon. G* is the free two-particle propagator. 

The nuclear structure effects are taken into account in Eq.(39) in terms of the nuclear charge 
radius rjv- This contribution has the form: 

^E h J r s SO p T (nS) = 2 < V«|AV£" • G ■ AV str \^ n >= ^^AE h F fs (nS)r 2 N G(0, 0). (40) 

The value of the reduced Coulomb Green function at zero arguments in the coordinate rep- 
resentation G(0, 0) is divergent. The reason of the appeared divergence lies in the expansion 
of the potentials in Eq.(40) at small relative momenta and further integration (40) at all 
values of relative momenta. To calculate the quantity G(0, 0) the dimensional regularization 
can be useful (29|, 30, 31]. Subtracting the iterative term 2 < ip n \ AV^ S ■ G? AV str \ipn > from 
Eq.(40) we obtain the following result: 
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One further contribution of the sixth order over a can be derived from the one-photon 
amplitude (Fig. 5 (a)) expanding the nuclear magnetic formfactor at small values of the rel- 
ative momenta. As a result the potential of the hyperfme interaction (3) in the coordinate 
representation gains the additional term 
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where tm is the nuclear magnetic radius. For the averaging the operator (43) over the 
Coulomb wave functions we use the following relation 
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and the value V 2 ^(0) = ^(0)/j, 2 (Za 
correction proportional to the nuclear magnetic radius: 



171 . |32| . As a result we obtain the important 
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TABLE I: Hyperfine structure of IS and 2S states in the ion of muonic helium (/i 2-He) + . ^12 = 
(8AE h f s (2S) - AE h f'(lS)). 



Contribution to HFS 


IS, meV 


2S, meV 


Interval A12 


Reference 


The Fermi energy 


-1370.725 


-171.341 





(6), [3,JJ 


Muon AMM correction of orders a 5 and a 6 


-1.598 


-0.200 





(7), [3,4 


Relativistic correction of order a 6 


-0.438 


-0.078 


-0.183 


(8), [3] 


One-loop VP contribution in I7 
interaction of order a 5 


-4.203 


-0.540 


-0.119 


(11)-(12) 


Two-loop VP contribution in I7 
interaction of order a 6 


-0.050 


-0.004 


0.016 


(14)-(16) 


One-loop muon VP contribution 
in I7 interaction of order a 6 


-0.052 


-0.007 





(11)-(12) 


One-loop VP contribution in the 
second order PT of order a 5 


-9.260 


-0.869 


2.305 


(23)-(24) 


Two-loop VP contribution in the 
second order PT of order a 6 


-0.105 


-0.010 


0.022 


(27)-(30) 


Nuclear structure correction of order a 5 


48.376 


6.047 





(33),(34) 


Contribution of VP and nuclear 
structure of order a 6 


0.760 


0.095 





(38) 


Nuclear structure correction 
of order a 6 


2.553 


0.272 


-0.377 


(41)-(42),(45) 


Nuclear structure and muon 
self-energy correction of order a 6 


-0.145 


-0.018 





(46), [2JJ 


Nuclear recoil correction 
of order a 6 


0.330 


0.038 


-0.026 


(47)-(48),[34j 


Summary contribution 


-1334.560 


-166.615 


1.638 





We have included in Table 1 the total nuclear structure contribution which is determined 
by expressions (41), (42) and (45) for IS and 2S energy levels at tm ~ r^. Let us write 
here also the corection of order a 6 connected with the nuclear structure and the muon self 
energy 



31 



^^str SE 



5 a(Za] 



71 



mxRz&E 



hfs 




0.145 meV 
0.018 meV' 



(46) 



where Rz is the Zemach radius. 

One part of the recoil corrections to the HFS is accounted in the calculation of the dia- 
grams in Figs. 3-4. Thus the leading order recoil contribution (ZoL){m\jm2) \w{m\lm-i)l\E } p S 
is contained in the potential (34). The recoil correction of order (Z a) 2 ( mi /m?)AEpf s for the 
ground state HFS in the hydrogen atom was calculated in Refs. [3l|, |33j and for the hyperfine 
interval A 12 in Ref.j34j. Using these results we can present the analytical expressions for the 
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G 



FIG. 5: Nuclear structure effects of order a in the one-photon interaction and the second order 
PT. G is the reduced Coulomb Green function. 



recoil corrections and their numerical values for the HFS of IS and 2S -states in the form 

■ 2 - 17 25 31C 



AE^ C S (1S) = {Zaf-^—/\E h F Ss {lS) 

777.1777.2 



12 + 3C + 72 + 



AE&(2S) = (Za) 



12 



VZa/ V 2 ' 2C 
2 - 265 821 809C 



0.330 meV, 



A££ /s (2S) 



777 1 777 2 



96 



+ 



+ 



+ hl2 M___,, +ln 



(za 



3 7 7C 

I h — 

2 2( 8 



96C 1152 
= 0.038 meV, 



(47) 



(48) 



where ( = 2m2^h/ r n p Z . 



V. SUMMARY AND CONCLUSION 

In this work various QED corrections, effects of the nuclear structure and recoil of orders 
a 5 and a 6 have been calculated for the hyperfine splittings of IS* and 25" energy levels in 
the ion of muonic helium (n\He) + . The investigation of the energy structure of IS* and 
2S states in this atom has the clear experimental prospect. Contrary to earlier performed 
studies of the energy spectra of light muonic atoms in Refs.0, 0, II] we use the three- 
dimensional quasipotential method for the description of the muon and helion bound state. 
All corrections considered here can be separated into two groups. The first group consists 
of the contributions which are specific for muonic helium ion. Primarily they are connected 
with the effects of the electron vacuum polarization. In our study these contributions are 
presented in the integral form and obtained numerically. The corrections known in the 
analytical form from the calculations of the hyperfine structure of muonium and hydrogen 
atom enter in the second group [3| • Numerical values of all corrections are presented in Table 
1. It contains also several basic references on the papers where the precision calculations of 
the HFS of simple atoms were considered. Other references can be founded in the review 
articles 0,3. 
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As mentioned above, the hyperfine structure of li ght exotic atoms was investigated on the 



basis of the Dirac equation many years ago in Refs. jlOi. Illj . The energies of the transitions 
(2S — 2P) were obtained for the muonic hydrogen and ion of muonic helium ([i\He) + . In this 
calculation only basic contributions to the HFS with the precision 0.1 meV were accounted. 
It follows from the Table 2 of Ref.[Hj], that the energies of the transitions Q~S\/i — z P\/2) and 
( 3 >S'i/2 — Z Pi/2) are equal correspondingly 1167.3 meV and 1334.1 meV resulting the definite 
value of the hyperfine splitting of 25-state: -166.8 meV. The total value of the HFS of 25* 
energy level entering in our Table 1 -166.615 meV is in good agreement with this result. So, 
the calculation of the HFS in the ion of muonic helium performed in this work improves the 
obtained earlier result in [ll[ for the 25* state by the calculation of order a 6 corrections and 
gives new result for the hyperfine splitting of IS state. The estimate of the next to considered 
order contribution over a has the form: a 3 ln(l/a)AEp fs (lS) « 0.003 meV. Despite the fact 
that all contributions in Table 1 are written with the accuracy 0.001 meV, the precision of our 
calculation of the HFS is not so high. The reason is that the nuclear structure correction of 
order a 5 has the largest theoretical uncertainty associated with the errors in the measurement 
of electromagnetic form factors for the nucleus \He. When the dipole parameterization for 
the form factors is used the value of the theoretical uncertainty is determined by the error 
of the nuclear charge radius: r^i^He) = 1.844 ± 0.045 fm. As a result the theoretical error 
can reach near ±1.5 meV for lS-level and ±0.20 meV for 2S'-level. Nuclear corrections to 
the HFS connected with the motion of the nucleons forming the nucleus 2 H, 3 H, 3 He of 



36]. Another source of the uncertainty is 



2$ [H, 0, H E IH, fl 13 • lt demands 



light hydrogen-like atoms were studied in Refs. 35, 
connected with the nuclear polarizability effect [37J, 
further investigation on the basis of the experimental data on the polarized lepton scattering 
by the nucleus \He. The value of the nuclear polarizability contribution to the HFS of 
the muonic helium ion can amount to several meV. The nuclear polarizability contribution 
should be considered in the combination with the nuclear corrections because both effects are 
associated with the interaction of the multinucleon system with electromagnetic field. The 
hyperfine splitting interval A12 has not uncertainties conditioned by the nuclear structure 
and polarizability corrections. So, the obtained value of the interval Aj.2 = 1.638 meV can 
be used for the check of QED predictions in the case of muonic helium ion with the accuracy 
0.01 meV. 
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